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This article examines the problem of periodic motions near collin- 
ear points of libration in the plane restricted problem of three bodies. 
With the aid of the Lyapunov methods, we demonstrate that around each 
collinear point of libration, a continuous family of periodic motions 
exists, dependent on one arbitrary constant. We examine the stability 
of the indicated periodic motions and show that they are unstable. 

Finally, we establish the orbital instability of the periodic motions 
under study. 

We shall proceed from differential equations, determining the 
movement of a particle (zero mass) near any collinear point of libration. 
To be definite, we shall examine point L 2 . The indicated equations can 
be written as follows: 


d 2 £ 

2n -^l = 

dn 

dt 2 

dt 

H 

d2 ? + 

2n-^- = 

d£2 

dt 2 

dt 

3q 


( 1 ) 


where £ and q --particle coordinates in a system of coordinates with its 
beginning at point L 2 , n 2 = 1 + jx and 12 is determined by formula 

ZQ = Pi 2 + + p ( pz2 + Jl) ’ 

where 

p ‘ 2 =( e + Q - iry)’ + v 

^ * ( e + ° + rfp :) 2 + ^ •) 

Finally, |a| is the distance of the point of libration to the center 
of gravity of two finite masses which are equal to 1 and |x, respectively. 
Since function!! is a holomorphic function of £ and r\ , at least in the 
direct proximity of the point of libration, it can be expanded into 
Taylor’s series, converging absolutely for all values of £ and rj , the 
moduli of which do not exceed known limits. This series can be 
written in the form 

n= k p ^ ? " i qr]Z ~ a ^ 3 + i a ^ 2 + _3 ^V 

(4) 

+ |( 3 t i 4 -yS 5 + 5 y£V - ^ + •••> 


(2) 


(3) 
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where 
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+ 4 f, 

q = - 1 -p. 

+ 2f, 

2f = 

1 

P ; 

1 

a = —5 — + 
P 10 

P 

- 3 I - 

Pio 

P 20 

P 10 


1 

M- 

1 

P 

p = 

"pT + 

P 5 20 ’ 

Y = ~b — + 
P 10 

O 

_ 

CL 



whereby 

Pio 


a 


P 

1 + P ’ 


P20 


a + 


1 

ttp * 


(5) 


( 6 ) 


If we retain only the terms of the second degree in our expansion of 
function fl, Equation (1) will be transformed into linear equations with 
constant coefficients. The determining equation of this linear system, 
as is well known, has the form 


x 4 - (p-q-4n 2 )x 2 -pq = 0. (7) 

This equation has the following roots: 

Xi, _ Xi, \j , _ \ u (i (8) 

where 

X 2 = 1 + n -f + V9f' -4f(l + (i), 

Xi = - 1 - M- + f + V9f 2 - 4f(l + n). 

\ Z always differs from zero and is positive. \i becomes zero only 
with (jl^O. Therefore, with all values of \x other than zero, Equation(7) 
possesses two purely imaginary and two real roots, one of which is 
positive and the other negative. The presence of one positive root of 
Equation (7) causes instability of collinear libration points in the 
Lyapunov sense, which is a well-known fact. 

Let us now examine the matter of the existence of periodic 
motions near collinear points. This problem is solved without 
difficulty, with the aid of Lyapunov's theorem alone, proof for which 
can be found in his well known essay 1 . This theorem can be for- 
mulated as follows: if the determining equation of a system of disturbed 

motion equations has two simple purely imaginary roots, whole 
multiples of which are not equation roots, and if the latter does not 
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possess zero roots, any time it is possible to find periodic series of 
a certain type which would formally satisfy the system of equations of 
disturbed motion, the latter, in fact, has a periodic solution presented 
by such series. 

Thus, to find a periodic solution, it is necessary to set up series 
of which all terms are periodic functions. Therefore, in the general 
case, it is impossible to establish the existence of periodic solutions 
a priori, for it is virtually impossible to determine all terms of an 
infinite series. 

Fortunately, there is at least one case whereby the existence of 
a periodic solution can be established a priori. This will be a case 
when equations of disturbed motion have a holomorphic integral, 
independent of time, of a definite structure. For example, if equations 
of disturbed motion have canonical form 


dx 3 9H dy 3 9H 

-dT = -977 -dT= tap <s=1 - 2 


• f n ) 9 


where H is a time -independent holomorphic function of variables x 3 and 
y 3 in which the terms of the lowest order generate a ground form, the 
equations have integral H= const. , and if the determining equation has 
a pair of simple purely imaginary roots \i and -XL, not possessing roots 
of the type mXi, where m=0, ± 1 , ±2, ± 3, . . . , the disturbed motion 
equations will, without doubt, have a periodic solution represented by 
periodic series of known form. 


Let us see what can be said of the existence of periodic motions 
near collinear libration points in the restricted problem of three bodies. 
As we have seen, the determining equation in this problem always 
has two simple purely imaginary roots and does not have (if p^O) zero 
roots or roots of type mXi, where m is a whole number. It is further 
known that system (1) has an integral independent of time. To become 
convinced that this integral has the required form, it is sufficient to 
reduce Equation (1) to canonical form, assuming 

> d dr) . 

* i= dF + n-e • 


Then Equation (!) will be replaced by the following system: 
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( 10 ) 


dt 


dt 


3H 

d £i 

3H 

~Hi ’ 

dt 

= ae ’ 

3H 

dr|i 

3H 

' 3t|i ’ 

dt 

3r| ’ 


where 


H=a- y (e 2 +t 1 2 )+n(eTn--nei)-y(ei 2 +tn 2 ). 


(ii) 


Since H is independent of time, Equation (10) has a holomorphic integral 
in which the terms of the lowest order assume ground form H 2 , of 
type 


H 2 = 2f£ 2 — fr| 2 + n(|r|i- r|^i) -y (^i 2 + r|i 2 ). 


( 12 ) 


Therefore, system (10) indubitably has a periodic solution, and 
consequently, such a solution will exist for system (1). 


Thus, we have established the existence of periodic motions 
near collinear points of libration, and it remains for us merely to 
find these motions. As Lyapunov 2 points out, for this we must proceed 
as follows. Assuming c to stand for arbitrary constant and T for 
series 


T = ^ (l + h 2 c 2 + h 3 c 3 


...) 


(13) 


with indefinite coefficients h 2 , h 3 , . . . , we shall insert in Equation (1) 
a new independent variable t to replace t by substituting 


t = 


2tt( t -t 0 ) 


(14) 


Constants h 2 , h 3 , . . . shall be selected in such a manner that transformed 
equation 


d 2 e 

- 2 nl. 

dr\ 

T 2 

30 

cTt 2 

2tt 

dr 

4tt 2 

ae ’ 

d 2 r] 

0 T 
-f 2n — 

2tt 

ii _ 

T 2 

30 

dT 2 

dT 

4tt 2 

an 


( 15 ) 
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will be satisfied by series 


£ = Xq = cx i + c Z x z + c 3 x 3 + . . . , } 

} (16) 
r| = y 0 = cyj + c 2 y 2 + c 3 y 3 + . . . , ) 

in which all coefficients, x 1} yi, x 2 , y 2 , ...» will be periodic functions 
of t with common period Zit . Since the existence of a periodic solution 
in our problem is known a priori, such series will indeed be found 
and will be absolutely convergent at all values of t and values of c 
satisfying condition 

| C | < C . (17) 

where c is a certain positive value differing from zero. Series (13) 
will also converge absolutely with |c| ^ C, and value T will be the 
period of the found motion (in respect to time). We have thus obtained 
a continuous series of periodic motions depending on two arbitrary 
constants, t 0 and c. Constant t 0 has no significant effect of the nature 
of the motion, since we can assume that the family of periodic motions 
depends solely on arbitrary constant c. 


Determining in practice functions xj, y lf x 2 , y 2 , . . . , we shall 
obtain for them expressions of form 


x i = Ai [cos r + sin t] , 

yi = Bx [cos r - sin tJ , 

x 2 = A 2 sin 2 t + C 2 , 
y 2 = B 2 cos 2t 

x 3 = A 3 [cos 3t - sin 3 t] , 

y 3 = B 3 [cos 3t + sin 3 tJ , 


( 18 ) 


where A x , B x , A 2 , B 2 , C 2 , 
F or example, 


Ai 


2n 

= ~ ’ 



are constants dependent only on \x . 


We shall not write out expressions for the remaining coefficients, in 
view of their awkwardness. 

Bet us now turn to an investigation of the stability of the periodic 
motions found by us. In the absolute sense, this problem is solved 
immediately and in the negative. Indeed, let us examine any of the 
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found periodic motions. It corresponds to a certain value of constant 
c, and consequently, a certain period T. Since the family of periodic 
solutions is continuous in respect to c, there exists an innumerable 
quantity of other periodic motions orbitally as close as desired to that 
under examination. The periods of these orbitally close motions will 
differ from the period of the motion in question, but they will differ 
in a degree as small as desired. Let us now picture two particles 
located at the initial moment in two periodic orbits as close as desired, 
and possessing initial velocity coordinate and component values as 
close as desired. Let us say that T and Ti are corresponding periods 
and let us assume, to be more definite, that T < Ti . It is obvious 
that after time interval T, the first particle will arrive at its initial 
state at that time when the second has not yet been able to reach that 
state. In other words, after time interval T, the second particle will 
lag somewhat behind the first and this gap will continue to grow with 
time. After a sufficiently large number of revolutions, our particles 
will be separated a distance comparable to the dimensions of their 
orbits and thus will cease being as close to one another as desired. 

This means that each periodic motion obtained by us is unstable in 
the Lyapunov sense. 

But instability of this type is not particularly interesting to us. 
Indeed, the behavior of the other motions in our problem is still an 
unknown, that is, motions not belonging to our periodic family. It 
may occur that trajectories of all remaining motions, as close as 
desired to periodic according to initial data, will always be as close 
as desired. This will mean that the periodic motion is orbitally 
stable. On the other hand, it may occur that motions will be found 
in our problem as close to periodic as desired according to initial 
data, the trajectories of which will not remain as close as desired to 
periodic orbit, which will signify that this periodic motion is also 
orbitally unstable. 

To elucidate this important matter of orbital stability of periodic 
motions near collinear points of libration we should examine Equation 

(15) , in which t is an independent variable. Indeed, in respect to t, 
all periodic motions of the family possess common period 2tt , as a 
consequence of which instability effect disappears (an effect occuring 
due to difference in periods of periodic orbits as close together as 
desired). We shall examine the problem of stability of any of solutions 

(16) of system (15). For this, we shall first insert new variables 
x and y in place of £ and rj, assuming 

t = x + x 0 , r| = y + y 0 . ( x 9) 
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The transformed equations will have the form 


d 2 x T 

dy 

T 2 

. dR 

dr^ 2n 2 tt 

dT 

4 t r 2 

' dx ’ 


dx 

rj-i2 

dR 

dT 2 2 IT 

dT “ 

4ir 2 

‘ 9y ’ 


where 

R(x, y) =Q(x + xo, y + y 0 ). 


( 20 ) 


( 21 ) 


Since all periodic orbits lie in an area where function Cl is 
holomorphic, function R(x, y) will also be a holomorphic function, at 
least for values of x and y,the moduli of which are sufficiently small. 
Therefore, function R(x,y) can be represented in the form of a series 
with coefficients periodic in relation to t. This series will be written 
as follows: 

R =-j P11 (t)x 2 + pi 2 (T)xy + ip 2 2 ( t )y 2 + . . . , (22) 


where the first coefficients will be determined by formulas 
Pii(t) = p-6axo + 12(3 xq-6Po + • • • » 


Pi 2 (t) = 3ay 0 -12(3x 0 y 0 + . . . , 

2 9 2 

P 22 ( T) = - q + 3ax 0 - 6(3x 0 + j (3y 0 + . . . 


(23) 


These coefficients are periodic functions of t, with common period 
2tt and, in addition, are dependent on parameter c. Since these 
coefficients are values of partial derivatives of holomorphic function 
fl(x + Xo, y+yo) with x=y=0, series (23) will be absolutely convergent 
in the area where this function is. holomorphic, and therefore 
coefficients Ph(t) , p 12 (t) , p Z 2(t), . . . can be represented in the form of 
series placed according to increasing degrees of parameter c and 
absolutely convergent at all values of t and with Icl ^ c. These 
series will have the form 
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Pii( T ) =p +cpi*j (t) + c 2 P ,2, (t) + . . . , 


P12(t) = cp*}^ (t) + cVl2 (T) + . . . , 

P 22 ( T) = - q + cpi ( T) + cp| ( T)+ . . . , 


k 1 ? ( r) 


( 24 ) 


where all coefficients pn( t), pi 2 (t), p 2 2 (t). . . are periodic functions of 
t with common period 2ir. Thus, our problem of periodic motion 
stability is reduced to a problem of stability of trivial equation 


x = 0, y = 0 

of a system with periodic coefficients (20). 


(ST) 


According to the general Lyapunov theory, we shall first examine 
the equations in variations corresponding to system (20). These 
equations in the variations will be written as follows: 


,2 

d x 
dT 2 


,2 

d y 




(25) 


^ y T dx T r / \ / \ 1 

^ + 2n i;7r [pu(r)x + p„(T )y ] 


and will represent a system of linear equations with periodic coefficients. 
The problem of stability (ST) of system (25) depends, as is known, on the 
nature of the roots of the characteristic equation of system (25). Since 
system (25) can be reduced to canonical form, its characteristic equation 
will have the form 3 

p 4 + Ap 3 + Bp 2 + Ap + 1 = 0, (26) 

and coefficients A and B will be holomorphic functions of parameter c 
and can be represented in the form of series located by increasing 
degrees of c and absolutely convergent at | c| ^ r 4 . 

To determine the nature of the roots of Equation (24), we shall 
see what we shall have with c = 0. Assuming in Equation (25) c = 0, 
we shall obtain a system with constant coefficients 
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( 27 ) 


d 2 x 

dT 2 


2n 

X” 


dy 

dT 


y ^ x - 0 

— ~w ~ u> 


d 2 y 

dT 2 


2n 

~ 


dx q 
dr + X 2 


-^- + tVy = o, 


the determining equation of which will have the form 

x 4 -— (p-q-4n 2 ) x 2 - — = 0. 

X 2 X 4 

Therefore, the roots of Equation (28) will be 



(28) 


(29) 


but system (27) can also be viewed as a system with periodic coefficients; 
therefore, for system (28), we can also build a characteristic equation 
of the same form as Equation (26). There is no need to set up the latter 
equation. The roots of the determining and characteristic equations 
of a system with constant coefficients are connected by formula 5 


P = 


e 21TX , 


(30) 


and since values of x are known, we immediately obtain corresponding 
values of p : 

2ttX i 2ttX i 

1, 1, e"* - , e X . (31) 

Hence, it is evident that with c = 0 Equation (26) possesses one root, 
the modulus of which is greater than one. Since A and B are holomorphic 
functions of c and since the above-mentioned root is simple, according 

to the theorem of the existence of implicit functions 6 , Equation (26) 

2TTX! 

with c X 0 will have a single root , tending to e“\ , when c tends to 

zero. Therefore, with sufficiently small values of |c| Equation (26) 
will obviously have a root the modulus of which is greater than one. 
According to Lyapunov's 7 results, we conclude that (ST) of system (25) 
and (ST) of system (20) will be unstable. 

This means that all periodic motions close to collinear points of 
libration in the restricted problem of three bodies also possess orbital 
instability, and the problem of stability of these periodic motions has 
thus been elucidated by us. 
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